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ABSTRACT 

A test of the equality of the parameters /\-j and A2 exponential 

random variables X-| and X2 is developed and evaluated for the case of a 
single sample from each population. The test is valid with or without 
time truncated testing although the test is less powerful under time 
truncated testing. The results are applicable to reliability studies. 
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I. INTRODUCTION 



A. PURPOSE 

The purpose of this thesis is to present a method for testing the 
equality of the scale parameters, and ^ 2 ’ O'*" exponential random 
variables X-j and X 2 whose respective density functions are 



v/hen time truncated testing of only one sample from each population is 
used. Time truncation consists of defining the new random variable Y 
where 



B. BACKGROUND 

Lentner and Buehler [Ref. 1] investigated the problem of making 
inferences about the relative magnitudes of the scale parameters of two 
untruncated exponential random variables X.| and X^ when sampling from 
each population was of size greater than one. In this case the sum of 
the random samples from each population is considered to be the observa- 
tion of a random variable Y.j whose distribution is gamma with parameters 

and n^. . The conditional distribution of the difference, Y-j -Y 2 
given the sum Y-j + Y 2 is dependent on the parameter Ai" A 2 • This result 
holds for any distribution from the exponential family, that is, a 
distribution of the form 



exp (- A^x.| ) X>0, x^ > 0 

exp (- X>0, X 2 > 0 




i = i ,2 



TyCy) = c(X)h(y) 
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This result was the basis for obtaining a confidence interval for the 



parameter A^ - A2‘ 

Lloyd and Lipow [Ref. 2] present methods for obtaining confidence 
intervals for the parameter \ in the case of a single population Y with 
an exponential distribution. These procedures utilize the chi-square 
distribution of the difference Y-.i-Y. i=l, ' • ’ ,n of the order statistics, 
0 =Yq<< Y^ ^ Y2^ Y3 ...•^Yp+1 for observations from an exponential 
di stribution. 

Cox [Ref. 3] uses the ranked observations from two exponential 
distributions and presents a modified F test for testing the hypothesis 
that both samples are from the same distribution, that is A-|=A2. The 
power, P [the statistic falls in rejection region], of this test was 
investigated by Gehan and Thomas [Ref. 4] for the case of time truncated 



toctinn nnri 1 1 .si7es 'JSinu L.uinuutt 



b i d i ! un . 






statistic here is t-j/t2 where 

t = j ^ij) - ( ^i + sum of observed values from sample i) 

^ (n^-v^) ( number of observed values from sample i ) 



v^. = number of truncated observations 
The test statistic has an approximate F distribution with degree of 
freedom 2(n-|-v.|) and 2(n2~V2). Their conclusion was that Cox's modified 
F test is nearly as powerful as the standard F test. The modified F test 
is not appropriate for the case of single samples since time truncation 
yields a zero divisor. 
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C. APPLICATIONS 



The results obtained herein are applicable to reliability studies of 
components assuming that the times to failure are distributed as exponen- 
tial random variables. The problem consists of choosing an alternate type 
of component to replace the type presently in use. In this case, high 
cost or large size might dictate the need to select an alternative with 
equal reliability when the items to be tested are in small supply and the 
time allocated is not sufficient for full testing. Another application 
for this procedure is the problem of making inferences when only one trial 
from each population was made, and outside considerations forced the 
termination of testing, for example, an electrical power failure during 
the initial testing of two limited production items. 
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II. DERIVATION OF PROCEDURE 



A. DISTRIBUTION OF THE SUM AND DIFFERENCE UNDER TIME TRUNCATION 
The random variable X with density function 
exp(-Ax) A>0,x^0 

is transformed to the random variable Y with density function 
exp(-Xy) 0<y<T 

exp(- XT) y=T 

when time truncated at time T. Using the simultaneous transformation 

Z,=V,+Y2 

The joint density function of and 1,^ may be found by use of Laplace 



T.rarir.Tnriiiarioiis. 



r i rs i I f ! fj. t f I s ! o y r I ' 1 1 u I • of V 



i ~ I 5 • I .1 IV.:** 



u (t) [^X exp(-Xt)J -u(t-T) 
+ u'(t-T) [^exp(-XT)] , 



Xex 



p(-Xt)] 



where u(t-a) is the unit step function at t=a, 



u(t-a)= 0 t<a 

-- 1 t>a 

I 

and u (t-a) is a unit mass function or Dirac delta function at t=a. 
u'(t-a) has the property 

fX 

I u ' (t-a )dt = 0 x<a 

-Co 

= 1 x>a 
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Now 1 (s^ , S2) = E l^exp (-s^Zi -S2Z2^ by definition and by 

independence of and this becomes 

(St ,s,) = f^ ■ 



'Z 7 '""I ’""2 

12 ' 



1 



For this case, using standard tables of Laplace transformations, 

f (s) = (A^+s)[] [j-exp(-(\^-+s )T)J + exp |^-(A^+s)tJ 

^i 

i = 1, 2. 

By multiplication of fv-iCsT+s^) and f^ (s -s^), then inverting the 

' 1 • ^ “2 1 ^ 

transformation, first with respect to s^ and then with respect to s^, 
the result, in terms of unit step functions and unit mass functions, is 
the joint density function of Z^ and Z^. This density is 

(1.1) (\^A2/2) exp[-(A^+ A 2 )z^/ 2 +(A 2 -Af;z^/ 2 ] 0<x.,<T 

1^2l<^l 



(1.2) A2®^P C(^2"V^'‘^2^i1 

(1.3) A^exp [(Ai-A2)T- A^zJ 

(1.4) exp [-(A^l+X^)!] 
v/hen A-|= A^= A, this becomes 



T<z^<2T 

IZ 2 U 2 T-Z 

T^z^<2T 

z =2T-z, 

2 1 

T^z^ < 2T 
Z2 =z,-2T 



z =2T, z =0 
‘ 2 
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(1.1.1) 


(X^/2) exp [-Az^ 


0<z.| <T, 
T^z^<2T, 


|Z2|<2T-Z.| 


(1.2.1.) 


exp [-Az^] 


Tiz^<2T, Z2=2 T-z^ 


(1.3.1.) 


exp [-z^X] 


TSz-]'‘^2T, z^-z^i -2T 


(•1.4.1.) 


exp [-2XT] 


CXI 

II 

N 


Z2=0 



To obtain the conditional density of given Z-j , the marginal density of 
Z^ is obtained by integration of the joint density vn'th respect to 
For example when T^z^<2T, f^ (z-j ) is 

^ ^ exp[-(\^+X2)z^/2+(\2-A^ ) z^/^ dz^ 

z -2T 

^ + A2GXP C(A2-A^)T-A2Z^1 + Aiexp[(A^-A2)T-A^z^^ . 

the marginal density of Z. is then 

I 

(1-5) [(A-|A 2 )/(A 2 "A] )] []exp(-A^z^ )-exp(-A 2 Z^ )} 0<z^<T 

[(AiA 2)/(A2-A^ )] [exp(-A2Z^ + X 2 T-A 1 T) 

-expC-A^z^-A^T+A T^j + X^exp [;(\2“A^ )T-A|Z^ 

+ A2exp[- A2Z^-(A^-A2 )t] 
exp[-(A^+ A^)t] 

and when A, this reduces to 

(1.5.1) (\^/2)[^exp (-Xz^)] (2z^) 

(1.6.1) (\2/2)[]exp (-Xz.,)] 2(2T-z^) 



( 1 . 6 ) 



(1.7) 



T<z^<2T 



Zi-2T 



0<z^<T 



(1.7.1) 



+2 Xexp(-Xz.| ) 
exp (-2 XT) 



T^z^<2T 

Zi=2T. 
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Finally, the density of given is obtained by division of (1.1) by 
(1.5), 

(1.8) 



( ^2" ^l1 


1 ( exp 


[(A 


2" >» ^)( 


:,)/2]\ 


V 2 j 


'\ 1- 


exp 




) 



0<z^<T 



and dividing (1.1) by (1.6) 



(1.9) 



^ 2 " 



exp 




^ Xi+X2 




A^exp J 


L"^ 2 ^r(^r^ 2 ^l] 


exp [ 


_-^i" 



-a2-A,)T] 



z,|<2T-z, 



Kz <2T 



The other portions of f^ ij will not be used in the development and are 
'..'hen A,=a^“ X , the cc.nditional density 



•Jt • 



(1.8.1) l/2z 



1 



iz kz, , 0<z^<T 



l/[2(2T-z^)+4/\] |Z2|<2T-Z^, T<z^<2T 

A heuristic ratio test, similar to a likelihood ratio test, ofX-| = ^ 
versus\,+A using the density function f is a test which calls for 

^ h ’^2 

rejecting the null hypothesis H^:A-| = A 2 when R defined as 

^^ 1 ^ 2 ‘^ 1^2 1 ^ 2 ^ 

^^ 1^2 



is less than or equal to a fixed number K and calls for accepting when 
R is greater than K. 
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In the problem being considered, if such a K exists then solving 



(X^/2)exp (-Xz-|) < K 

(\^A 2 / 2 ) exp [(A 2 -A^)z^/ 2 -(X^+A 2 )z^/ 2 ] 

for ^2“^1 ^ region of the form [(z,.2j,): lz2-z,|<g(k.z,)} 

called the critical region. Any observation falling in this region 
calls for rejecting the null hypothesis. The probability of an observa- 
tion falling in the critical region is called the power of the test. 

If time truncated testing had not been used and a trial 
solution of jz^l^ (l-ot)z^ for the critical region is assumed, the power of 
the test is 



C( 



■ oo 



f -7 7 dz^dz, + 



-d-cK)zi 



0 / - 7 - 






/I ^1^2 ^ ^ 

I 

and if A-|=A 2 =X, this reduces to 

^0 [‘^^ll ^ • 

When X-|'"A 2 , the probability that an observation falls in the critical 
region is, with truncation, 

T fz fT r-l(l-«c)z 

0 ] (Ljz/V^^ V ] 0 ] -z, 




r2T/(2-o()r: 

)t J 



f 2T/(2-0()('2T-z 

(l-c<)z J T Jz,-2T 



r 2T/(2-<x)r-(l-<x)z, 

1 f. . dz.dz, + 1 I fz,z/^2'^l 



r2T/(2-(x) r 
)t =2T-z, 



f^ dz^dZi 

-2 = 2^-1 

This is 1 minus the value listed in Table 1 



13 




i 




The use of time truncation maps the region 



A = 
onto 
B 









= -^^2! 



2T-Z1 , Tsz,s|L 



IJ , (1-<x)zi < [Z 2 |< 2T-Z, , T^z, < 2l 7. 
^ 2-x ) 



Clearly B is contained in A and A is contained in the critical region, thus 
under the null hypothesis the probability of observing [z^j ^(l-(x)zi is 
less than or equal to oc. This is the probability of rejecting the null 
hypothesis when is true. Figure 1 displays a lower bound for the 
power of this test whena =.3 and T= 21 . 

If testing is not time truncated, a power function which is monotoni- 
cally increasing in the absolute value of (A2"Ai )Z-| is available. The 
test statistic is Zp given Z., , and the decision rule is: reject A-|=A2 

if Z2^(l-(x)Zi. 

The power of the test in this case is 



r (i-x)z, 

1- ) ^z Iz 

^ -(l-oc)z-, ^2'^1 2 

_ i_ exp [(X2 "Ai ) (^1 "^^1 "^1 )/2] “6xp{^“(A2"Ai ) (z-j -c<z.j+z.| )/2j 

l-exp[-(\2-A^ )z^] 

sinh [(A2-A1) (l-o<)zi/2] 

^ sinh [(A 2 -A 1 ) z-|/2] 

whenA-| = A2 the power is o< as may be verefied by L. ‘ Hospital's rule. 

Figure 2 displays the power function without time truncation whencx=.l and 
for selected values of Zy 
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The probability of type II error is the probability that an observation 
falls in the compliment of the critical region, |z 2 l <(1 -tx)zi J , 

when Hq is not true. This probability is 



f 2 T/( 2 -c<) r 
^ 0 ^ 



2 T/( 2 -c<) r(l--(x)Zi 

-(l-ix)z.| 1' 2 



rOD r z,-(a2T)/(2-<x) 

+ \ \ fyfy dZ dZ 

^2T/(2-oc) Ja2T/(2-(x)-z.| Wz ^ ' 

(1.10) = 4(1-(x)(AiX2)/ (2-(x)Xi+ c<A2_ Jz-dc)A2+c<X^J 

^ exp[-XiT+X20<T/(2-o<)] ^ exp[-X2T+AixT/(2-(x)] ^ ^ 

(2-o<)X^-HxA 2 '^2^^ "T^A 2 -^X^ 



Table 1 shows this probability for selected values of X-| >X 2 > T andcx. When 
T becomes arbitrarily large, (1.10) reduces to 4(l-(x) (XiX-?)/ r( 2 -o()X,+cxX 2 T 

[(2-(x)X2+X-|0c] . 



B, SUMfvlARY OF TEST PROCEDURE 

The procedure for testing the hypothesis X-|=X 2 against the hypothesis 
^A'^l“^2 significance cx 

1 . Observe Y-| and 

2. Form Z^=Y^+Y 2 and Z 2 =Y^-Y 2 

3. Compute (l-o<)Z.j 

4. If ^ 2 "^ (1"C<)Z.| accept . 

If Z 2 ^(l-oc)Z.| reject Hq . 
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Table I. Probability of Type II Error 
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Figure 1. Power Function Lover Bound With Time Truncation. 
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Figure 2. Pov/er Functici Without Time Truncation 



III. ESTIMATION OF A2" 



A. POINT ESTIMATION OF X^ 

In section II, the density of I 2 given was derived for the case 
where neither X.jnor X^are time truncated. This is the case 0<Z^<T and 
the density is given by equation 1.8. The expected value of given Z^ 

7 ( ^ exp[(X 2 -A^)(z 2 -z^)/z] 

1\ 2 / l-exp[-(\2-Xi )zil 

2 l+exp[-(A2-A])Zil 2_ 

' l-exp[-(A2-A,)z,] ’^2-Al 

-1 

I / \ \ \ . in\ n t f\ \ \ 

z-j cotii Mrt2"Ai;^I/^ ' VA2 A|/- 

From this an estimate of the parameter X 2 "A| obtained by solving 

2^=2^ coth 1 ^X 2 " A] )zi/^ - 2 /(X 2 "Ai) 

forA 2 "A]* Numerical computations using Newton's method are appropriate 
with this equation. A typical example might be z^=2, Z2=-1.32. Then 

-1.32=2 coth (X 2 "A] ) -2/(X2"Ai) 

and the value of A 2 “A-| which satisfies this is X2"Ai“"3. 

This result is also valid for the case T^Z-|<2T and Z^?^ 2T-Z^ since 
neither X^nor X^have been time truncated in this case. 

For the case in which X^or X^are time truncated, a point estimate of 
X 2 ~\-] is not available. But a confidence interval forA 2 "Xi niay be obtained. 



is then 
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B. CONFIDENCE INTERVALS FOR \2~h 

A 100(1-0()% confidence interval for a parameter is an interval on the 
real line, (a,b), such that the estimator of that parameter has a l-c< 
probability of being observed in the interval (a, b) and the unknov/n 

parameter will lie in the interval. 

For the exponential random variable X with scale parameter X, the 
unbiased estimator for Xis the observation X. A 100(l-(x)% confidence 
interval for Xis 



( -ln{l-0(/2) , 


-In (a/2) \ 


V X 


X ) 



To obtain a (l-(x) 100% confidence interval for X2"^l neither X.jnor 



confidence intervals for X-| and ^ 2 * Then by the transiormation X ^ 2 "^! » ^ 
confidence interval with confidence greater than or equal (l-<x) 100% is 

- 1n(l-^/2) + lnV2 X 



When either X^or X^are time truncated, a (l-o<) 100% confidence interva 



again using the transformation X=\ 2 -.\i » assuming and Y 2 =T, a 

confidence interval with confidence greater than or equal (l-(x) 100% is 

( ln(l-^/2) , - In^ + ln^/2 \ 

\ y T y / 



X 2 is time 



truncated, first obtain separately (1-3) 100%, 3 -l-A-^x, 




for X is 
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C. APPLICATIONS OF CONFIDENCE INTERVALS 

An ad hoc procedure for accepting the hypothesis Xi=A 2 based on 
confidence intervals for X 2 "X] would be to compute the confidence interval 
based on the observed values of the random variables Y-j and '< 2 - Then if the 
point 0 is in the interval obtained, accept \i=X 2 - This procedure 
corresponds favorably with the procedure proposed in Section II. 

To illustrate this method and its correspondence with the proposed 
procedure of Section II, consider the observations 

Y^=T=1 and Y 2 = .6666 

Here Z^=1.6666 and Z2=.3333. Whan(X=.3, (1 -oc)Z^=l .1 666 which is greater 
than Z 2 , hence accept ^^q-X-| = X 2 * 

Pop }^=1-A^=.165, an 83.5% confidence interval for A-, is (0, -In. 165) = 
(0, 1.800) and an 83.5% confidence interval for X 2 is 

( (-In. 9175)/. 666, (-In .0825)/. 666) ) 

= ( .129, 3.76 ) . 

Now, a confidence interval for X 2 "X-| with confidence greater than or equal 
to 70% is 

( -1.671 , 3.76 ) 

which includes the point zero, so the hypothesis X-| = X 2 should be accepted. 
This agrees with the test result from Section II. 
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IV. APRIORI KNOWLEDGE OF A 



A. CHOICE OF TEST TIME T. 

When the value of the scale parameter is known apriori, the experimenter 
may reduce the probability of type II error by a more advantageous choice 
of the test time T. If is considered to be the standard, hence known 
apriori, the probability of type II error, given by equation 1.10, can be 
bounded. Specifically, the terms involving A 2 T niay be bounded by choosing 
T appropriately. For example, if the constraint 

exp [-X 2 T] 

is to be imposed, then T must be greater than -In-^/ • 

For reliability studies, the scale parameters are usually known from 
previous data and hence can be of use for this type or test. 
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V. AREAS FOR FURTHER STUDY 



A. EXTENSION OF RESULT TO SAMPLES GREATER THAN ONE 

In order to extend the results obtained thus far to tests with more 
than one sample point from each population, it would be necessary to 
determine the joint density of the sum of the Z-j^- random variables and the 
sum of the l 2 ^ random variables. Here the independence of random variables 
within populations and the independence between populations of the random 
variables makes the computations easy since 



m ^ 

E z ^ % z 

i li > i 2i 1 ’ ^2 



(s) 



where an equal number, m, of observations are taken from each population. 

A». -J i-i 4* -I ir. w *I T nrr 1;.-,. | r.-F vFcin H ". n n Y'Pi'Z':\T hP TP. 

^ - ^ - 

test 

X! Z > (l-<x)X z 
i i 

and reject the hypothesis H^= X-| = !^2 inequality was observed. Here 

one needs the density function of this test. Without truncation, the 
density of the sum E Z^i given the sum X Z.| . , is the density whose 

Laplace transform is 



exp(-5 E Z„ 

which by independence is 



m 

7T E 

l=i 







exp(-sZ2i) 
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When produce of m Laplace transformations of uniform 

(-z^^. , ) density functions, hence it is 

exp(sz^^. ) - exp(-sz^.) 

- 




i Tf 

sm i=i 



sinh(^^1 i ) 



-li 



For m larger than 2, this is difficult to invert in order to obtain the 
density function. 

Another possibility is to consider the index sets 



J = fi 


: Z <t| 




li ^ 


K= (i 


: TZ^.^2T , 


r 


1 


L= 


^2i""'"“-|i J 


M = 


: Z2.=Z,.-2t} 


N = ^i 


: Z,.=2t] 



< 2T-Z 



I 

lii 



and then consider a likelihood ratio test 



^ ^2i ’^7i ^^2i^li I 

Using the independence of the paired observations this ratio is the product 
over the respective index sets of the densities for the region corresponding 
to the index sets. Solving this for z^^ yields a critical region of the 
form 



iej 



>..+X: 

i6K 



Z^..+Z_i Z +c^ 
2 1 2i 1 
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This suggests using the previously mentioned test 

S a-- (!-«) 2 z,i 

i i 

but with the modification of introducing scaling factors for pairs belonging 
to index sets L and M. Such a modification might be investigated as part of 
a future study. 



25 



VI. CONCLUSIONS 



The test prescribed in this thesis has significance level less than 
or equal 0 (, if time truncation is used, and significance level 0 (, if time 
truncation is not used, when testing the equality of parameters for random 
variables with exponential distributions. As the parameter X2~ A.| moves 
away from the value zero, the power function increases. 

When time truncation is not used, the power function is a monotonical ly 
increasing function of the absolute value of {X2"Ai)^1* 

The test prescribed is applicable when testing with only a single 
sample from each population and it can be used instead of Cox's modified 
F test in the case of single samples when time truncation is used. 

Thp Mcp nf rnnfidence intervals described in this thesis is an 
extension of the procedure of Lentner and Buehler to the case of single 
samples from each population. This procedure may be used in the case of 
time truncated testing; the procedure of Lentner and Buehler may not. 
However, the confidence intervals used here are not exact, while those 
of Lentner and Buehler are. 
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